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Abstract 

We study three-flavor neutrino oscillations in the early universe in the presence of 
neutrino chemical potentials. We take into account all effects from the background 
medium, i.e. collisional damping, the refractive effects from charged leptons, and 
in particular neutrino self- interactions that synchronize the neutrino oscillations. 
We find that effective flavor equilibrium between all active neutrino species is es- 
tablished well before the big-bang nucleosynthesis (BBN) epoch if the neutrino 
oscillation parameters are in the range indicated by the atmospheric neutrino data 
and by the large mixing angle (LMA) MSW solution of the solar neutrino problem. 
For the other solutions of the solar neutrino problem, partial flavor equilibrium may 
be achieved if the angle #13 is close to the experimental limit tan 2 #13 < 0.065. In 
the LMA case, the BBN limit on the v e degeneracy parameter, \^ v \ < 0.07, now ap- 
plies to all flavors. Therefore, a putative extra cosmic radiation contribution from 
degenerate neutrinos is limited to such low values that it is neither observable in the 
large-scale structure of the universe nor in the anisotropies of the cosmic microwave 
background radiation. Existing limits and possible future measurements, for ex- 
ample in KATRIN, of the absolute neutrino mass scale will provide unambiguous 
information on the cosmic neutrino mass density, essentially free of the uncertainty 
of the neutrino chemical potentials. 
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1 Introduction 



The cosmic matter and radiation inventory is known with ever increasing 
precision, but many important questions remain open. The cosmic neutrino 
background is a generic feature of the standard hot big bang model, and its 
presence is indirectly established by the accurate agreement between the cal- 
culated and observed primordial light-element abundances. However, the exact 
neutrino number density is not known as it depends on the unknown chemical 
potentials for the three flavors. (In addition there could be a population of 
sterile neutrinos, a hypothesis that we will not discuss here.) The standard 
assumption is that the asymmetry between neutrinos and anti-neutrinos is of 
order the baryon asymmetry rj = (% — ng)/n 7 ~ 6 x lCT 10 . This would be 
the case, for example, if B — L = where B and L are the cosmic baryon and 
lepton asymmetries, respectively. While B — L = is motivated by scenarios 
where the baryon asymmetry is obtained via leptogenesis [1], there are models 
for producing large L and small B [2-5]. 

In order to quantify a putative neutrino asymmetry we assume that well before 
thermal neutrino decoupling a given flavor is characterized by a Fermi-Dirac 
distribution with a chemical potential \i v , f u (p, T) = [exp (E p /T — £„) + 1] , 
where £„ = \i v jT is the degeneracy parameter and E p ~ p since we may neglect 
small neutrino mass effects on the distribution function. For anti-neutrinos the 
distribution function is given by £ p = 

A neutrino chemical potential modifies the outcome of primordial nucleosyn- 
thesis in two different ways [6]. The first effect appears only in the electron 
sector because electron neutrinos participate in the beta processes which de- 
termine the primordial neutron-to-proton ratio so that n/p oc exp(— £ e )-Q 
Therefore, a positive £ e decreases Y p , the primordial 4 He mass fraction, while 
a negative £ e increases it, leading to an allowed range 



compatible with £ e = (see Refs. [7-10] and Sec. 5). A second effect is an 
increase of the neutrino energy density for any non-zero £ which in turn in- 
creases the expansion rate of the universe, thus enhancing Y p . This applies to 
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all flavors so that the effect of chemical potentials in the or v T sector can be 
compensated by a positive £ e . Altogether the big-bang nucleosynthesis (BBN) 
limits on the neutrino chemical potentials are thus not very restrictive. 

Another consequence of the extra radiation density in degenerate neutrinos 
is that it postpones the epoch of matter-radiation equality. In the cosmic 
microwave background radiation (CMBR) it boosts the amplitude of the first 
acoustic peak of the angular power spectrum and shifts all peaks to smaller 
scales. Moreover, the power spectrum of density fluctuations on small scales 
is suppressed [11,12], leading to observable effects in the cosmic large-scale 
structure (LSS). 

A recent analysis of the combined effect of a non-zero neutrino asymmetry on 
BBN and CMBR/LSS yields the allowed regions [13] 

- 0.0K £ e < 0.22, |^ lT |<2.6, (2) 

in agreement with similar bounds in [14,15]. These limits allow for a very 
significant radiation contribution of degenerate neutrinos, leading many au- 
thors to discuss the implications of a large neutrino asymmetry in different 
physical situations. These include the explanation of the former discrepancy 
between the BBN and CMBR results on the baryon asymmetry [16] or the ori- 
gin of the cosmic rays with energies in excess of the Greisen-Zatsepin-Kuzmin 
cutoff [17]. In addition, if the present relic neutrino background is strongly de- 
generate, it would enhance the contribution of massive neutrinos to the total 
energy density [18,19] and affect the flavor oscillations of the high-energy neu- 
trinos [20] which are thought to be produced in the astrophysical accelerators 
of high-energy cosmic rays. 

The limits in Eq. (2) ignore neutrino flavor oscillations, an assumption which 
is no longer justified in view of the experimental signatures for neutrino oscil- 
lations by solar and atmospheric neutrinos. For zero initial neutrino chemical 
potentials, the flavor neutrinos have the same spectra so that oscillations pro- 
duce no effect. This is true up to a small spectral distortion caused by the 
heating of neutrinos from e + e~ annihilations, an effect which is different for 
electron and muon/tau neutrinos and which causes a small relative change in 
the final production of 4 He of order 10 -3 [21]. This relative change is slightly 
enhanced by neutrino flavor oscillations [22,23]. In the presence of neutrino 
asymmetries, flavor oscillations equalize the neutrino chemical potentials if 
there is enough time for this relaxation process to be effective [24]. If flavor 
equilibrium is reached before BBN, then the restrictive limits on £ e in Eq. (1) 
will apply to all flavors, in turn implying that the cosmic neutrino radiation 
density is close to its standard value. As a consequence, it is no longer neces- 
sary to use the neutrino radiation density as a fit parameter for CMBR/LSS 
analyses, unless one considers exotic models with decaying massive particles. 
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The effects of flavor oscillations on possible neutrino degeneracies have been 
considered in [20] , where it was concluded that flavor equilibrium was achieved 
before the BBN epoch if the solar neutrino problem was explained by the large- 
mixing angle (LMA) solution. The LMA solution is favored by the current solar 
neutrino data. Thus, it was concluded that in the LMA case a large cosmic 
neutrino degeneracy was no longer allowed. 

We revisit this problem because the flavor evolution of the neutrino ensemble 
is more subtle than previously envisaged if medium effects are systematically 
included. Contrary to the treatment of Ref. [20], the refractive effect of charged 
leptons can not be ignored, and actually is one of the dominant effects. While 
the background neutrinos produce an even larger refractive term, its effect 
is to synchronize the neutrino oscillations [25,26] which remain sensitive to 
the charged-lepton contribution. Still, equilibrium is essentially, but not com- 
pletely, achieved in the LMA case so that our final conclusion is qualitatively 
similar to that of Ref. [20] . 

One counter-intuitive subtlety is that the neutrino self-potential actually can 
suppress oscillations in a situation where the excess of neutrinos in one flavor 
is exactly matched by an excess of anti-neutrinos in another flavor. In this 
case the synchronized oscillation frequency is zero so that oscillations begin 
only once the cosmic expansion has diluted the self-term. Therefore, one can 
construct cases where flavor equilibrium is not achieved before BBN even in the 
LMA case. However, this is only possible for specially chosen initial conditions 
where is equal for all flavors, but the absolute signs may be different. For 
the purpose of deriving limits on however, this case is equivalent to the 
one where equilibrium is achieved, the only important point being that |^| is 
approximately equal for all flavors at the BBN epoch. 

Another subtlety appears in a full three-flavor analysis. We show that achiev- 
ing equilibrium in the LMA case does not depend on the value of the mixing 
angle #13, to which strict limits from reactor experiments apply. Moreover, for 
the non-LMA solutions of the solar neutrino problem, partial flavor equilib- 
rium may be reached if the angle #13 is small but non-zero. 

In Sec. 2 we set up the formalism to study primordial flavor oscillations. We 
then turn in Sec. 3 to the primordial flavor evolution of a simplified system 
where v e mixes maximally with one other flavor. This two-flavor case will il- 
lustrate many of the important subtleties of our problem. Then we turn in 
Sec. 4 to realistic three-flavor situations which involve yet further complica- 
tions. In Sec. 5 we finally derive new limits on the degeneracy parameters and 
summarize our findings. 
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2 Neutrino flavor oscillations in the early universe 



In order to study neutrino oscillations in the early universe we characterize 
the neutrino ensemble in the usual way by generalized occupation numbers, 
i.e. by 3x3 density matrices for neutrinos and anti- neutrinos as described in 
[27,28]. The form of the density matrices for a mode with momentum p is 



Pip, t) 



( Pee Pen Per \ 

Pfie PfMfM PflT 
\ Pre Pt/j, Ptt j 



piP, t) 



^ Pee Pfie Pre ^ 
Pen Pun Pi~n 

\ Per PfiT Ptt j 



(3) 



where overbarred quantities refer to anti-neutrinos. The diagonal elements 
of the density matrices correspond to the usual occupation numbers of the 
different flavors. The definition of the density matrix for anti-neutrinos is 
"transposed" relative to that for neutrinos, allowing one to write the equations 
of motion in a compact form [27]. 

The equations of motion for the density matrices relevant for our situation of 
interest are [27,29] 




3m? 



E + P - P , P P 



E + p- p) ,p p 



+ C[ Pp ] , 

+ C\p p \ , (4) 



where we use the notation p p = p(p,t) and [•, •] is the usual commutator. 
Further, M 2 is the mass-squared matrix in the flavor basis; in the mass basis 
it would be diag(raf , ra|, m§). The diagonal matrix E represents the energy 
densities of charged leptons. For example, E ee is the energy density of elec- 
trons plus that of positrons. The density matrix p is the integrated neutrino 
density matrix so that, for example, p ee is the v e number density while p ee 
is the v e number density. The term proportional to (p — p) is non-linear in 
the neutrino density matrix and represents self- interactions. Finally, C[-] is the 
collision term which is proportional to G|. We have neglected a refractive term 
proportional to the neutrino energy density which is much smaller than the 
(p — p) term in our present situation where the neutrino asymmetries are as- 
sumed to be large. On the other hand, we have neglected the usual term which 
is proportional to the charged-lepton asymmetries. This asymmetric term is 
always negligible; at early times (high temperatures) it is negligible compared 
to the E term, while at temperatures near n/p freeze out (T ~ 1 MeV) it 
is negligible compared to the vacuum term M 2 /2p if the mass-squared dif- 
ferences coincide with those characterizing the atmospheric neutrino and the 
solar neutrino LMA oscillations. 
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In the expanding universe we need to substitute d t — > d t — Hpd p with H the 
cosmic expansion parameter. We have taken this into account rewriting the 
equations in terms of comoving variables, as listed in the Appendix. We solve 
these equations numerically, calculating the evolution of p and p on a grid of 
neutrino momenta. 



3 Two-flavor oscillations 



3. 1 Equations of motion for the polarization vectors 



In order to develop a first understanding of flavor oscillation in the early 
universe it will be useful to study first a two-flavor situation involving u e and 
Vp. The usual relation between flavor and mass eigenstates is 



v e = cos 9 v\ + sin 9 v 2 , 
Up = — sin + cos 9 v 2 . 

The 2x2 density matrices are 



(5) 



p(p, t) = p*) = \ [ p °o°> *) + CT • p 0°' *)] , 

p(p,t) = (f 1 

\ Pefj, rmi ■ 



P (p,t) + a-P(p,t) 



(6) 



Here, Oi are the Pauli matrices while P(p,t) and P(p,t) are the usual polar- 
ization vectors for the neutrino and anti-neutrino modes p, respectively. We 
normalize p and p to a Fermi-Dirac distribution with zero chemical poten- 
tial, f(p) = [exp (p/T) + so that for instance f Va (p,t) = paa(p,t)f(p)- 
For the polarization vectors, the equations of motion are given by the usual 
spin-precession formula 



d t P P = + 



d t P P = - 



Am 2 D 8V2G fPj7 ; 
B — — — E ee z 



2p 

Am 2 
~^p~ 



B - 



3m w 



\y/2G F p 
3771^ 



X P r 



+ v / 2G F (P-P) x P P + C[P P ] , 



X P r 



+ v / 2G F (P-P)xP p + C[P p ] , (7) 



where Am 2 = m 2 , — m 2 and B = (sin 29, 0, — cos 29) with the vacuum mixing 
angle. Further, z is a unit vector in the z-direction, E ee the electron-positron 
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energy density, and P and P are the integrated polarization vectors. 

We will see that a detailed treatment of the collision terms is not crucial. 
Therefore, we approximate them with a simple damping prescription of the 
form 

C[PJ = -A.PJ (8) 

where the transverse part Pj consists of the x-y-projection of P p . We use 
the damping functions as, for instance, in [28,30], but with slightly modified 
coefficients. We neglect a small dependence on £, i.e. we use the same damping 
functions for neutrinos and antineutrinos. We have checked that our results are 
insensitive to the inclusion of the repopulation function, i.e. Pq and P as given 
for instance in [28,30]. Moreover, we did not include the neutrino heating from 
e + e~ annihilations, since it is small even in the presence of degeneracies [31]. 

3.2 Evolution of simple flavor asymmetry 

As a first case we consider a simple situation where initially electron neutrinos 
do not have a chemical potential (£ e = 0) while muon neutrinos are asymmetric 
(^n = —0.05). We use a rather small asymmetry so that the anticipated 
equilibrium state £ e = ^ will be close to the BBN limit on £ e in Eq. (2). The 
quantitative evolution with a large £ M is the same. Moreover, we use maximal 
mixing and Am 2 = 4.5 x 10~ 5 eV 2 as suggested by the LMA solution of the 
solar neutrino problem. 

In order to illustrate the relative importance of the various contributions to the 
equation of motion we first calculate the primordial evolution of the integrated 
polarization vectors P and P for pure vacuum oscillations, including only the 
Am 2 /2p terms in Eq. (7), without any medium effects whatsoever. Our initial 
conditions £ e = and £ M = —0.05 imply that there are more v e than i.e. P 
points initially in the +z direction, while there are more than v e , i.e. P 
points initially in the — z direction. In the upper panel of Fig. 1 we show the 
evolution of P z as a function of cosmic temperature; the evolution of P z is 
similar, except with a negative initial value. The oscillations begin when the 
expansion rate has become slow enough at about T ~ 30 MeV. The oscillation 
frequencies are different for different modes, leading to quick decoherence and 
thus to an incoherent equal mix of both flavors. Evidently flavor equilibrium 
is reached long before n/p freeze-out at T ~ 1 MeV. 

As a next step we include the neutrino self-potential proportional to P — P. 
The evolution is shown as a dotted line in the upper panel of Fig. 1. It is 
now dominated by the effect of synchronized oscillations, i.e. the self-potential 
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forces all neutrino modes to follow the same oscillation [26] . Because P and P 
point in opposite directions, the common or synchronized oscillation frequency 
^synch is very much smaller than a typical Am 2 /2p. Explicitly we have [26] 

= jif / M-tr 1 ■ (Pp + Pp) s w (9) 



where I = P — P and I is a unit vector in the direction of I. Numerically our 
initial conditions imply p cf[ ~ 132 T. 

Next we add to the vacuum term the effective potential caused by the e ± back- 
ground [see Eq. (4)], without including the self-term. The e ± potential always 
points in the z-direction, thereby suppressing flavor oscillations in the usual 
way. This effect is shown in the lower panel of Fig. 1 (solid line) . P stays frozen 
at its initial value up to T ~ 3 MeV where the medium potential becomes 
unimportant compared to the vacuum term. At these temperatures collisions 
are also about to become unimportant. Therefore it is not surprising that 
including collisions (dotted line) does not dramatically change the evolution 
of P z . 

However, this impression is misleading as can be seen from Fig. 2 where we 
show the evolution of P x for the same cases as in the lower panel of Fig. I. 
Without collisions (solid line) the evolution of P is a simple turning from the 
z- to the ^-direction. The potential adiabatically disappears, leading to the 
usual MSW-type evolution. In this case the final result corresponds to equal 
numbers of both flavors, yet a coherent flavor superposition. With collisional 
damping (dotted line) the re-component is damped, leaving little flavor co- 
herence of the final state. Therefore, with or without collisional damping we 
reach "flavor equilibrium" in the sense of equal densities of both flavors, but 
only collisions ensure the damping of the otherwise large transverse part of 
the final P. 

Finally we study how oscillations evolve in the presence of all effects: vacuum, 
e ± background, neutrino self-potential and collisions. The results are shown in 
Fig. 3 for different strengths of the self-potential. In one run it was switched 
off ("No Self"), it was taken at full strength ("Self"), or it was suppressed by a 
factor 10 4 or 10 5 . When the self-interaction term is close to its real value (solid 
line), again all modes oscillate in synch so that we see a combination of the 
synchronized oscillations and the MSW-like effect of the background medium. 
However, collisions ensure that no large transverse component survives — see 
the corresponding evolution of P x and P y in Fig. 4. Almost perfect flavor 
equilibrium is achieved with or without self-interactions. 

This conclusion is not to be taken for granted. When the self-interaction term 
is artificially adjusted to be of comparable strength to the other terms (vacuum 
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or e ± background), the evolution shows complicated features (Fig. 3). After 
some initial conversion, P z remains constant for a long time, only to reach 
equilibrium at a much later time. This behavior is not easy to understand 
as it depends on all the ingredients of our calculation, i.e. the background 
potential, the neutrino self-potential and damping by collisions. We believe 
that the long phase where the polarization vector essentially stands still is 
caused by the synchronized oscillation frequency becoming very small. This is 
strictly the case in a situation where the chemical potentials for two flavors 
are equal but opposite — see Sec. 3.3 below. We believe that in the present 
case a combination of decoherence of some of the modes and collisions drives 
the system into a state where the synchronized oscillations frequency vanishes 
so that oscillations stop until the cosmic expansion has diluted the neutrino 
self-term enough to eliminate this effect. 



3.3 Equal but opposite asymmetries 

Our discussion thus far suggests that for the parameters of interest it is generic 
to achieve flavor equilibrium before the BBN epoch, even though the loss of 
coherence may not always be complete. However, this conclusion depends on 
the simple initial conditions that we have used thus far. The different neutrino 
flavors may show large asymmetries, yet the total cosmic lepton asymmetry 
could still be small if initially £ e = — £ M , which corresponds to conservation of 
the lepton number L e —L^. These initial conditions imply that P = — P so that 
the self-potential term (P — P) is large while for each mode (P p + P p ) = [26]. 
This implies, in turn, that the synchronized oscillation frequency u> sy nch — 0. 
Therefore, in this case the synchronization effect caused by the self-potential 
prevents flavor oscillations entirely, at least until the self-term becomes weak, 
long after the BBN epoch. While flavor equilibrium is not achieved in this 
case, from the start we have |£ e | = so that the BBN limits on |£ e | would 
apply to both flavors. 



4 Three-flavor oscillations 



4-1 Oscillation parameters 

We now turn to the more interesting case of three-flavor neutrino oscillations. 
The neutrino flavor eigenstates u e , u^, and v T are related to the mass eigen- 
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states via the mixing matrix 
/ 



C12C13 S12C13 S13 

-S12C23 - C12S23S13 C12C23 - S12S23S13 S23C13 
•S12S23 — C12C23S13 — C 12 S 2 3 — Sl2 c 23 s 13 C23C13 



(10) 



Here = cos#y and s^- = sin^- for ij = 12, 23, or 13, and we have assumed 
CP conservation. The set of oscillation parameters is now five- dimensional (see 
for instance [32]), 



Am sun 


= Am 2 21 




= a^1 2 


^sun 


= #12 


#atm 


= #23 


#13 





1 1 

m 3 ~ m 2 



;n) 



We do not perform a global analysis of all possible values of these parameters, 
but fix them to be in the regions that solve the atmospheric and solar neutrino 
problems [32,33]. In particular we take Am 2 tm = 3 x 10 -3 eV 2 and maximal 
mixing for 8 atm from the former, while from the solar analyses we consider the 
following values for Am 2 un in eV 2 : 4.5 x 10~ 5 , 7 x 10" 6 , 1 x 10~ 7 , 8 x 10~ n 
for the Large Mixing Angle (LMA), Small Mixing Angle (SMA), LOW and 
Vacuum regions, respectively. For the angle # sun we take the approximation of 
maximal mixing for all cases except SMA where we use 9 sun = 1.5°. 



4-2 Simple three flavor case 



We begin with the simplified situation where initially only the muon neutrinos 
are asymmetric (£ M = —0.1). We perform a full three-flavor calculation, but 
for now set #13 = 0. The evolution of the neutrino asymmetries is shown 
in Figs. 5 and 6 for the LMA and LOW cases, respectively, both with and 
without the neutrino self-interactions. For this choice of oscillation parameters 
the three-flavor oscillations effectively separate as two two-flavor problems for 
the atmospheric and solar parameters, respectively. The oscillations caused by 
the largest Am 2 are effective at T ~ 20 MeV, as soon as the ^ background 
disappears completely. The presence of the self-term causes only a slight delay 
in the equilibration of £ M and £ r . 

The oscillations due to Am 2 un and # sun are effective only when the vacuum term 
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overcomes the e ± potential. In the LMA case, the conversions takes place above 
T ~ 1 MeV, leading to nearly complete flavor equilibrium before the onset of 
BBN. For the LOW parameters the synchronized oscillations just start at that 
epoch. The presence of the neutrino self-potential does not significantly change 
the picture in the LMA case while for the LOW case one clearly observes the 
phenomenon of synchronized oscillations. For the SMA and Vacuum regions 
primordial oscillations involving u e are not effective before BBN if 9 13 = 0. 

4-3 Non-zero 13 mixing 

The angle 13 is restricted to the approximate region tan 2 9 13 < 0.065 (see 
for instance [34]) from a combined analysis of solar, atmospheric and reactor 
(CHOOZ) data. However, a small but non-zero #13 does modify the oscillation 
behavior. This effect is shown in Fig. 7 for the LMA case and different values 
of #13. Even small values of #13 lead to conversion to the electronic flavor at 
larger temperatures and enhance flavor equilibration if 9 sun is in the LMA 
region. 

The effect of a non-zero 9 13 is more important if the solar neutrino problem 
is solved by oscillations with parameters in a region other than LMA. This 
can be seen in Fig. 8, where we have calculated the evolution in the limit 
Am 2 un ~ and 9 sun ~ 0, thus corresponding approximately to the SMA, 
LOW and Vacuum regions. For values of tan 2 # 13 close to the limit discussed 
at the beginning of this section, neutrino oscillations lead to almost complete 
flavor equilibration, while for smaller #13 angles the conversion is only partial. 

4-4 Can the self-term prevent flavor equilibrium? 

In the two-flavor case of Sec. 3.3 we saw that for equal but opposite asym- 
metries of the two flavor distributions the synchronized oscillation frequency 
was strictly zero, suppressing oscillations entirely before the BBN epoch. In 
the three flavor case the situation is more complicated so that again we need 
to raise the question if there is a special configuration of initial neutrino dis- 
tributions which suppresses flavor oscillations by the neutrino self-potential. 

In analogy to Sec. 3.3 we first consider situations where the asymmetries of 
two flavors are equal but opposite, while the third asymmetry is strictly zero. 
In this case one expects that oscillations between the asymmetric flavors are 
suppressed by the self-potential while oscillations into the third flavor are 
unimpeded. We show the numerical results for two specific cases. In Fig. 9 we 
take £ M = — £ r = —0.1 and £ e = 0. The "atmospheric" oscillations between 
and v T are indeed blocked in the presence of self-interactions while the "solar 
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oscillations", here the LMA case and tan 2 #13 = 0, proceed at the appropriate 
temperature. In Fig. 10 we take £ r = and £ e = — £ M = 0.1. Both the atmo- 
spheric and solar oscillations take place and equilibrate the flavors essentially 
as in the previous cases with a simple initial condition. 

The question remains if one could devise special initial conditions such that 
\C, T \ and are large while |£ e | is small, and that this system avoids flavor 
oscillations by reducing the synchronized oscillation frequency to a very small 
value. In this case our new limits would not apply. 

We do not believe that such a case can be constructed. To argue in favor 
of this claim we first note that the most general initial condition consists of 
thermal equilibrium distributions which are diagonal in flavor space, i.e. in 
the weak- interaction basis the initial 3x3 density matrices for neutrinos and 
anti-neutrinos are diagonal and characterized by Fermi-Dirac distributions. 
All off-diagonal elements would be quickly damped by reactions which involve 
muons and electrons and thus "measure" the flavor content of the participating 
neutrinos. Therefore, the most natural initial condition for our problem is 
indeed characterized by £ e , £ M and £ T . 

The general equations of motion Eq. (4) reveal that the quantity which de- 
velops in a synchronized fashion is p — p, in agreement with the discussion 
of 2x2 synchronized oscillations in Ref. [26] where the picture of polarization 
vectors and the associated "internal magnetic fields" of the system was used. 
Therefore, the synchronized equation of motion is obtained by subtracting the 
two lines of Eq. (4) and integrating over all modes, 



where we have dropped the collision term. The last term, of course, is iden- 
tically zero. The second term vanishes also as long as the density matrices 
are diagonal in the flavor basis since the matrix E is diagonal in that basis. 
Only the first term, involving the mass matrix, represents a "force" such as 
to move p — p away from being diagonal in the flavor basis, i.e. which causes 
synchronized flavor oscillations. This force is identically zero if the matrix 



is proportional to the unit matrix. As the matrices p p and p p are initially 




+ V2G F [(p-p),(p-p)] , 



(12) 




(13) 
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diagonal and given by Fermi-Dirac distributions, the integral can be solved 
explicitly, leading to an expression proportional to 









\ 


(' ' ] 




i 


2 








M 



(14) 



If all £ are initially equal, then the neutrinos are already in flavor equilibrium 
and indeed oscillations trivially do not operate. 

The only nontrivial configuration where synchronized oscillations proceed with 
a vanishing frequency is when one of the chemical potentials is opposite in sign 
to the others. We have checked numerically that indeed no flavor conversion 
arises in this case as long as the self-potential is large, i.e. as long as we are 
in the synchronized regime. In the two-flavor case this corresponds to the 
situation discussed in Sec. 3.3. 

As in the two-flavor case we conclude that it is possible to avoid flavor equi- 
librium by specially chosen initial conditions, but these conditions require 
|£e| = = |£t|- Therefore, the strict BBN limits on |£ e | apply to all flavors. 



5 New limits on neutrino degeneracy 



We conclude that in the LMA region the neutrino flavors essentially equili- 
brate long before n/p freeze out, even when # 13 is vanishingly small. For the 
other cases the outcome depends on the magnitude of #13. In the LMA case it 
is thus justified to derive new limits on the cosmic neutrino degeneracy param- 
eters under the assumption that all three neutrino flavors are characterized by 
a single degeneracy parameter, independently of the primordial initial condi- 
tions. We do not derive the corresponding limits for the other solar neutrino 
solutions, since they would strongly depend on the value of a non-zero #13. 
However, if that angle is close to the experimental limit, the bounds that we 
describe would approximately apply. 

We first note that the energy density in one species of neutrinos and anti- 
neutrinos with degeneracy parameter £ is 



7tt 2 



" 120 




(15) 



It is clear that the BBN limit will imply £ 1 for all flavors so that the mod- 
ified energy density and the resulting change of the primordial helium abun- 
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dance Y p will be negligibly small. If there are additional relativistic species, 
such as sterile neutrinos or majorons, then Eq. (2) will simply apply to all the 
active neutrinos 



Therefore, the only remaining BBN effect is the shift of the beta equilibrium by 
£ e . We recall that Y p is essentially given by n/p at the weak- interaction freeze- 
out, and that n/p oc exp(— £ e ) ~ 1 — £ e where the latter expansion applies 
for |£ e | < 1. Therefore, AY P ~ -F p (l - Y p /2% ~ -0.21 £ e . Modifications of 
Y p by new physics are frequently expressed in terms of the equivalent number 
of neutrino flavors AN V which would cause the same modification due to the 
changed expansion rate at BBN. If the radiation density at BBN is expressed in 
terms of N u , the helium yield can be expressed by the empirical formula AY P = 
0.012 AN U [35]. Therefore, the effect of a small f e on the helium abundance is 
equivalent to AN V ~ — 18£ e . A conservative standard limits holds that BBN 
implies |AiV„| < 1 which thus translates into |£ e | < 0.057. 

A more detailed recent analysis reveals that the measured primordial helium 
abundance implies a 95% CL range N v = 2.5±0.8 or AN V = -0.5±0.8 [8,13]. 
We conclude that the BBN-favored range for the electron neutrino degeneracy 
parameter is at 95% CL 



If all degeneracy parameters are the same, then this range applies to all flavors. 

It should be noted that the actual limit we obtain on the neutrino degeneracy 
depends on the adopted BBN analysis. For instance AN U could be as high 
as 1.2 when the primordial abundance of lithium is used instead of that of 
deuterium [36]. At any rate, a limit of |£ e | < 0.1 seems rather conservative 
and does not modify our conclusions. 

Using |£| < 0.07 as a limit on the one degeneracy parameter for all flavors, 
the extra radiation density is limited by (Ap vD )/ ' p vV < 3 x 0.0021 = 0.0064, 
i.e. AN U < 0.0064. If the same radiation density were to be produced by the 
asymmetry of one single species, this would correspond to |£| < 0.12. 

For comparison with the future satellite experiments MAP and PLANCK that 
will measure the CMBR anisotropies, it was calculated that they optimistically 
will be sensitive to a single-species £ above 0.5 and 0.25, respectively [37]. 
However, with proper consideration of the degeneracy with the matter density, 
and the spectral index, n, a more realistic sensitivity is £ ~ 2.4 and 0.73, 
respectively [38]. Turning this around we conclude that our new limits are 
so restrictive that the CMBR is certain to remain unaffected by neutrino 



If I <0.22. 



(16) 



f e = 0.03 ±0.04 . 



(17) 
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degeneracy effects so that |£| can be safely neglected as a fit parameter in 
future analyses. 

If our new limits apply the number density of relic neutrinos is very close to its 
standard value. Therefore, existing limits and possible future measurements 
of the absolute neutrino mass scale, for example in the forthcoming tritium 
decay experiment KATRIN [39], will provide unambiguous information on the 
cosmic mass density in neutrinos, free of the uncertainty of neutrino chemical 
potentials. 

Note added: Very recently the authors of refs. [40,41] have also analyzed the 
equilibration of neutrino asymmetry from flavor oscillations, providing further 
analytical insight and confirming our conclusions. 
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A Evolution equations 

In this appendix we list in detail the evolution equations for the neutrino and 
anti- neutrino density matrices in the two-flavor case as in Eq. (5), while the 
generalization to the three-flavor case that we consider in Sec. 4 is straight- 
forward. 

In our treatment of primordial neutrino oscillations we use the following di- 
mensionless expansion rate and momenta 

x = mR , y = pR , (A.l) 

where R is the universe scale factor and m an arbitrary mass scale that we 
choose to be 1 MeV. The neutrino and anti-neutrino density matrices are 

p(x,y) = 1 -iP (x,y) + *.P(x,y)]= 1 -( K ^+^ ( A - 2 ) 
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with the derivatives 



i \ tt dp , . dp 



(A.4) 



The initial conditions are for large temperatures as follows. For density ma- 
trices normalized to f cq (y) = (e y + l) -1 and initial degeneracies £ a = — and 
£a = — £0 (for flavor neutrinos z/ a and i/g) 



Po(y) = 
^(1/) = 
Po(y) = 



feq(y,£a) + fcqiV,^) 

/eq(y) 
/eq(l/, ~£q) + /cqtj/, 
/eq(?/) 



/cq (?/) 

where f cq {y,0 = l/[exp(y - + !]■ Finally, 

^(y) = n(j/)=^x(s/) = P„(2/) = o 



(A.5) 



(A.6) 



The components of the neutrino polarization vectors evolve as 



dPo 

dx 
dP* 
dx 

dPy 

dx 



dP^ 
dx 

dP 
dx 

dP x 
dx 



(x,y) = R a (x,y) + Rp(x,y) , 

^ y)= {0Tx cos26 -^x) P ^ y) - Wx P ^ y) ' 

(Am 2 V \ D 
cos26> I P x (x,y) PJx,y) 
2pHx Ex) K ,yJ Ex yK ,UJ 

Am 2 



>,y) 



2pEx 
Am 2 



sin26 l P z (x,y) , 



sin 26 P y (x, y) + R a (x, y) - Rp (x, y ) , 



2pEx 

(x,y) = R a (x,y) + Rp(x,y) , 

( Am 2 Vi \ — D — 

^y) = -{ m - x ^ 2d -jrx) p ^ y) ~ Wx p ^ y) ' 
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dx 



dP 2 



( Am 



\2pHx 
Am 2 



2 Vi \ D 



Hx 



dx 



t ------ sin 20 P z (x,y) , 

sin 26 P y (x, y) + R a (x, y) - Rp(x, y) 



2pHx 
Am 2 
2pHx 



where 



(A.7) 



R a (x,y) = 2— 



(A.8) 



with R a fi given by the same expressions with £ — > — £. However, our numeri- 
cal results do not change significantly if the total neutrino and anti-neutrino 
number densities are taken to be constant (dP /dx = dP /dx = 0). 

The different terms in Eqs. (A.7)are given as follows. The vacuum oscillation 
terms are proportional to 



Am 2 _ 10 10 M P (Am 2 \ J_ x"_ 

2 p Hx ~ 2^7/i V ev 2 ) VP y 



(A.9) 



where M P = 1.221, p = (x/m) 4 p tot , and p tot is the total energy density of the 
universe. The / + /~ background with I — /i,e the charged leptons (for z/ r -f M 
and v e -v^ T oscillations, respectively) is described by 



Hx 



8^2 10 5 M P G F 1 y _ 

~Fk 3d KP1++P1-) 



3 J87T/3 m 2 w VP x 4 



(A.10) 



where Gf = 1.1664 and my/ = 80.42. Collisions are described by damping 
terms where, according to our calculations, 



D~2 x (4 sin 4 9 W - 2 sin 2 9 W + 2)F , 
D~2 x (2sin 4 ^ + 6)F , 



(A.ll) 



for v t -Vij, and v e -v^ T oscillations, respectively, where all collision terms are 
proportional to 



Fh_ = M P G 2 F C(3) y mcd _L_ j/_ 
^ 3 7T 3 ycV3 \/P xi 



(A.12) 
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and C(3) ~ 1.20206 and y mcd ~ 3.15137. 

Finally one has to add the contribution from the neutrino-antineutrino back- 
ground [27-29] 



^ (J(x)-J(s))-^ (v(x)+V(xj) 
(j( x )_J( x )) + ^E (U(x)+U(a; 



x P(x,j/) , 

xP(x,y), (A. 13) 



iJx ^ ' Hx 

for dP/dx and dP/dx, respectively. Here we have defined 

J(z) = ^J duu2 P (^ M ) > = ^J duu3 P ( X ' M ) ' ( A - 14 ) 

and the same for J and U. Moreover, 
Ksym _ V2 10 11 M P G F 1 



^ ^8^/3 >/P ^ 2 

cos 2 ^ — A . (A. 15) 



y sym _ 8^ 10 5 MpG F 2n 1 y 



Hx 
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Pure vacuum 
Vacuum + Self 
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Vacuum + e - potential 
-0.002 -Vacuum + e ± potential + Collisions 



10 1 

T (MeV) 

Fig. 1. Evolution of the z-component of the integrated polarization vector P for 
different situations as explained in the text. 
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Fig. 2. Evolution of the x-component of the integrated polarization vector P for the 
same cases as in the lower panel of Fig. 1. 
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Fig. 3. Evolution of the z-component of the integrated polarization vector P in the 
presence of the vacuum term, the e + e~ potential, collisional damping, and varying 
strength of the neutrino self-interaction. 
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Fig. 4. Evolution of the x and y components of the integrated polarization vector P 
in the presence of all effects and full-strength neutrino self-interactions, correspond- 
ing to the solid line in Fig. 3. 
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Fig. 5. Evolution of the neutrino degeneracy parameters for the LMA case, #13 = 0, 
and the initial values £ e = £ T = and £^ = —0.1. 
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Fig. 6. Evolution of the neutrino degeneracy parameters for the LOW case, #13 = 0, 
and the initial values £ e = £ T = and ^ = —0.1. 
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Fig. 7. Evolution of the neutrino degeneracies for different values of tan 2 #13 in the 
LMA case for the initial conditions £ e = £ r = and £ M = —0.1. 
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Fig. 8. Evolution of the neutrino degeneracies for different values of tan 2 613 in the 
limit Am 2 un ~ and 6 sun ~ for the initial conditions £ e = £ r = and £ M = —0.1. 
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Fig. 9. Evolution of the neutrino degeneracies if initially £ e = and £ r = — = 0.1, 
with or without neutrino self-interactions. 
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Fig. 10. Evolution of the neutrino degeneracies for the LMA case with #13 = 
with initial conditions £ r = and £ e = — £ M = 0.1, with or without neutrino 
self-interactions. 
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